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Abstract 

Ultrahigh energy (UHE) extragalactic protons propagating through cosmic mi- 
crowave radiation (CMB) acquire the spectrum features in the form of the dip, 
bump (pile-up protons) and the Greisen-Zatsepin-Kuzmin (GZK) cutoff. We have 
performed the analysis of these features in terms of the modification factor. This 
analysis is weakly model-dependent, especially in case of the dip. The energy shape 
of the dip is confirmed by the Akeno-AGASA data with = 19.06 for d.o.f. = 17 
with two free parameters used for comparison. The agreement with HiRes data is 
also very good. This is the strong evidence that UHE cosmic rays observed at ener- 
gies 1 X 10^^ eV - 4 X 10^^ eV are extragalactic protons propagating through CMB. 
The dip is also present in case of diffusive propagation in magnetic field. 
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1 Introduction 



The nature of signal carriers of UHECR is not yet established. The most 
natural primary particles are extragalactic protons. Due to interaction with 
the CMB radiation the UHE protons from extragalactic sources are predicted 
to have a sharp steepening of energy spectrum, so called GZK cutoff [1]. 

There are two other signatures of extragalactic protons in the spectrum: dip 
and bump [2,3,4,5]. The dip is produced due to p + 7cmb — ^ p + e"*" + e~ 
interaction at energy centered by ~ 8 x 10^*^ eV. The bump is produced by 
pile-up protons which loose energy in the GZK cutoff. As was demonstrated in 
[3], see also [5], the bump is clearly seen from a single source at large redshift 
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but it practically disappears in the diffuse spectrum, because individual 
peaks are located at different energies. 

As will be discussed in tfiis paper, the dip is a rehable feature in the UHE 
proton spectrum (see also [6] - [8]). Being relatively faint feature, it is however 
clearly seen in the spectra observed by AGASA, Fly's Eye, HiRes and Yakutsk 
arrays (see [9] and [10] for the data). We argue here that it can be considered 
as the confirmed signature of interaction of extragalactic UHE protons with 
CMB. 
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Fig. 1. The HiRes data [11] on mass composition. The measured atmospheric depths 
of EAS maximum Xmax at E > I x 10^^ eV (triangles) are in a good agreement 
with QGSJet-Corsika prediction for protons. 

The measurement of the atmospheric height of EAS maximum, Xmax, in the 
HiRes experiment (see Fig.l) gives another evidence of the proton composition 
of UHECR at > 1 X 10^^ eV. Yakutsk [12] and HiRes-Mia [13] data also 
favour the proton composition at > 1 x 10^* eV, though some other methods 
of mass measurements indicate the mixed chemical composition [14]. 

At what energy the extragalactic component sets in? 

According to the KASCADE data [15], the spectrum of galactic protons has 
a steepening at £^ ~ 2.5 x 10^^ eV (the first knee), helium nuclei - at £^ ~ 
6 X 10^^ eV, and carbon nuclei - at ~ 1.5 x 10^^ eV. It confirms the rigidity- 
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dependent confinement with critical rigidity Rc = Ec/Z fa 3 x 10^^ eV. Then 
CR galactic iron nuclei are expected to have the critical energy of confinement 
at Ec ^ I X 10^^ eV, and extragalactic protons can naturally dominate at 
£■ > 1 X 10^^ eV. This energy is close to the position of the second knee (Akeno 
- 6 X 10^^ eV, Fly's Eye - 4 x 10^^ eV, HiRes - 7 x 10^^ eV and Yakutsk - 
8 X 10^^ eV). The detailed analysis of transition from galactic to extragalactic 
component of cosmic rays is given in [16]. It favours the transition at E ^ 

1 X 10^^ eV. The model of galactic cosmic rays developed by Biermann et al 
[17] also predicts the second knee as the "end" of galactic cosmic rays due to 
rigidity bending in wind-shell around SN, produced by Wolf-Rayet stars. The 
extragalactic component becomes the dominant one at energy ~ 1 x 10^^ eV 
(see Fig.l in [17]). 

Below we shall analyze the features in UHE proton spectrum using basically 
two assumptions: the uniform distribution of the sources in the universe and 
the power-law generation spectrum. We shall discuss how large-scale inho- 
mogeneities in source distribution affect the shape of the features. We do not 
consider the possible speculations, such as cosmological evolution of sources. In 
contrast to our earlier works [6,7,16], we do not use here the model-dependent 
complex generation spectrum: the modification factor method allows us to use 
more general power- law spectrum with an arbitrary 7^. 

2 Bump in the diffuse spectrum 

The analysis of the bump and dip is convenient to perform in terms of modi- 
fication factor [3]. 

The modification factor is defined as a ratio of the spectrum Jp{E), with all 
energy losses taken into account, to unmodified spectrum J^"™, where only 
adiabatic energy losses (red shift) are included, 

r,iE) = J,{E)/j;^-{E). (1) 

For the power-law generation spectrum oc Eg'^" from the sources without 
cosmological evolution one obtains the unmodified spectrum as 

jn^) = ^(7. - moE-'^ J dz-{i + zro^\ (2) 



where the observed energy E and emissivity Cq are measured in GeV and 
GeV/Mpc^yr, respectively. The connection between d,t and dz is given by 
usual cosmological expression (see e.g. [7]). The fiux Jp{E) is calculated as in 
[6] with all energy losses included. 
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Fig. 2. Modification factor for the power-law generation spectra with 7^ in the 
range 2.0 - 2.7. Curve rj = 1 corresponds to adiabatic energy losses, curves rjee - to 
adiabatic and pair production energy losses and curves r/tot - to all energy losses. 

In Fig. 2 the modification factor is shown as a function of energy for two 
spectrum indices 7^ = 2.0 and 7^ = 2.7. They do not differ much from each 
other because both numerator and denominator in Eq. (1) include factor E~'^^ . 
Let us discuss first the bump. We see no indication of the bump in Fig. 2 at 
merging of rjee{E) and r]tot{E) curves, where it should be located. The absence 
of the bump in the diffuse spectrum can be easily understood. The bumps 
are clearly seen in the spectra of the single remote sources [3]. These bumps, 
located at different energies, produce a fiat feature, when they are summed up 
in the diffuse spectrum. This effect can be illustrated by Fig. 5 from Ref. [3]. 
The diffuse flux there is calculated in the model where sources are distributed 
uniformly in the sphere of radius -Rmax (or -Zmax) • When Zmax are small (between 
0.01 and 0.1) the bumps are seen in the diffuse spectra. When radius of the 
sphere becomes larger, the bumps merge producing the fiat feature in the 
spectrum. If the diffuse spectrum is plotted as E^Jp{E) this flat feature looks 
like a pseudo-bump. 



3 Dip as a signature of the proton interaction with CMB. 

The dip is more reliable signature of interaction of protons with CMB than 
GZK feature. The shape of the GZK feature is strongly model-dependent: 
it is more fiat in case of local overdensity of the sources, and more steep in 
case of their local deficit. It depends also on fluctuations in the distances 
between sources inside the GZK sphere and on fluctuations of luminosities 
of the sources there. The shape of the dip is fixed and has a specific form 
which is difficult to imitate by other mechanisms. The protons in the dip are 
collected from the large volume with the radius about 1000 Mpc and therefore 
the assumption of uniform distribution of sources within this volume is well 
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Fig. 3. Predicted dip in comparison with the Akeno-AGASA data. 

justified. In contrast to this well predicted and specifically shaped feature, 
the cutoff, if discovered, can be produced as the acceleration cutoff. Since the 
shape of both GZK cutoff and acceleration cutoff is model-dependent, it will 
be difficult to argue in favour of any of them. The problem of identification of 
the dip depends on the accuracy of observational data, which should confirm 
the specific (and well predicted) shape of this feature. Do the present data 
have the needed accuracy? 

The comparison of the calculated modification factor with that obtained from 
the Akeno-AGASA data, using ■jg = 2.7, is given in Fig. 3. It shows the 
excellent agreement between predicted and observed modification factors for 
the dip. 

In Fig. 3 one observes that at E < 1 x 10^^ eV the agreement between calcu- 
lated and observed modification factors becomes worse and at E < 4x 10^^ eV 
the observational modification factor becomes larger than 1. Since by defini- 
tion ri{E) < 1, it signals about appearance of another component of cosmic 
rays, which is most probably galactic cosmic rays. The condition t] > 1 means 
the dominance of the new (galactic) component, the transition occurs at higher 
energy. To calculate for the confirmation of the dip by Akeno-AGASA data, 
we choose the energy interval between 1 x 10^^ eV (which is somewhat arbi- 
trary in our analysis) and 4 x 10^^ eV (the energy of intersection of r)ee{E) and 
rjtotiE)). In calculations we used the Gaussian statistics for low-energy bins, 
and the Poisson statistics for the high energy bins of AGASA. It results in 
= 19.06. The number of Akeno-AGASA bins is 19. We use in calculations 
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Fig. 4. Predicted dip in comparison with the HiRes data. 

two free parameters: 7^ and the total normalization of spectrum. In effect, 
the confirmation of the dip is characterised by = 19.06 for d.o.f.=17, or 
X^/d.o.f.=1.12, very close to the ideal value 1.0 for the Poisson statistics. 

In Fig. 4 the comparison of modification factor with the HiRes data is shown. 
The agreement is also very good: = 19.5 for d.o.f. = 19 for the Poisson 
statistics. 

The good agreement of the shape of the dip riee{E) with observations is a 
strong evidence for extragalactic protons interacting with CMB. This evidence 
is confirmed by the HiRes data on the mass composition (see Fig. 1). 

The dip is also present in case of diffusive propagation in magnetic field [18]. 

4 Extragalactic iron nuclei as UHECR primaries 

Does modification factor for iron nuclei differ from the proton dip? 

We calculated the modification factor for iron nuclei, assuming that that Fe 
nuclei are the heaviest ones accelerated in the sources, and considering the 
propagation of Fe nuclei with energy losses taken into account. The resulting 
flux is given only for primary iron nuclei, without secondary nuclei produced 
during propagation (more details will be presented in [19]). 

The energy losses for Fe are dominated by adiabatic energy losses up to 4.5 x 
10^^ eV, from where on e"'"e~ energy losses dominate. For energies E > 1.7 x 
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Fig. 5. Modification factor for iron nuclei in comparison witli that for protons. 
Curve rj = 1 corresponds to adiabatic energy losses. Proton modification factors are 
given given by curve 1 (adiabatic and pair production energy losses) and by curve 2 
(total energy losses). Iron modification factors are given by curve 3 (adiabatic and 
pair production energy losses) and by curve 4 (with photodissociation included). 

10^° eV photodissociation becomes the main source of energy losses [21,19]. 
According to this dependence the energy spectrum of iron nuclei is oc E~'^3 up 
to ii^i ~ 1 X 10^^ eV, where the first steepening begins. The second steepening, 
caused by iron-nuclei destruction, occurs at energy i?2 ~ 1 x 10^° eV. For 
lighter nuclei the steepening (cutoff) starts at lower energies. Therefore the 
cutoff of the nuclei spectra occurs approximately at the same energy as the 
GZK cutoff, though the physical reason for these two cutoffs is different: while 
the latter (GZK) is caused by starting of photopion production, the former 
(nuclei) - by transition from adiabatic to pair-production energy losses [20]. 

In Fig. 5 the modification factors for iron nuclei are shown as function of energy 
in comparison with modification factors for protons. Comparison with Fig. 3 
clearly shows that even small admixture of iron nuclei in the primary extra- 
galactic flux upsets the good agreement of the proton dip with observational 
data. 

5 Discussion and conclusions 

There are three signatures of UHE protons propagating through CMB: GZK 
cutoff, bump and dip. 

The energy shape of the GZK feature is model dependent. The local excess 
of sources makes it flatter, and the deficit - steeper. The shape is affected 
by fluctuations of source luminosities and distances between the sources. The 
cutoff, if discovered, can be produced as the acceleration cutoff (steepening 
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below the maximum energy of acceleration). Since the shape of both, the GZK 
cutoff and acceleration cutoff, is model-dependent, it will be difficult to argue 
in favour of any of them, in case a cutoff is discovered. 



The bump is produced by pile-up protons, which are loosing energy in pho- 
topion interactions and are accumulated at low energy, where the photopion 
energy losses become equal to that due to pair-production. Such bump is dis- 
tinctly seen in calculation of spectrum from a single remote source. In the 
diffuse spectrum, since the individual peaks located at different energies, a 
flat spectrum feature is produced. 

The dip is the most remarkable feature of interaction with CMB. The protons 
in this energy region are collected from the distances ~ 1000 Mpc, with each 
radial interval dr providing the equal flux. All density irregularities and all 
fluctuations are averaged at this distance, and assumption of uniform distri- 
bution of sources with average distances between sources and average lumi- 
nosities becomes quite reliable. The dip is confirmed by Akeno-AGASA and 
HiRes data with the great accuracy (see Figs 3 and 4). As one can see from 
Fig. 5, presence of even small fraction of extragalactic heavy nuclei in the 
primary flux upsets this agreement. 

We interpret the excellent agreement of the calculated dip with the observa- 
tions as an independent evidence that observed primaries at energy 1 x IQ-*^^ — 

4 X 10^^ eV are extragalactic protons. This evidence is the complementary one 
to the direct measurements (now contradictive) of chemical composition. 

At energy ii^ < 4 x 10^'^ eV the modification factor from Akeno data exceeds 1, 
and it signals about dominance of another cosmic ray component, most prob- 
ably the galactic one. It agrees with transition from galactic to extragalactic 
component at the second knee ~ 1 x 10^^ eV. This conclusion is confirmed 
by the recent HiRes data on mass composition (see Fig. 1) and indirectly by 
the KASCADE data (see [16] for the detailed analysis). 

Are there alternative explanations of the dip? The conservative one (sec e.g. 
[22]) is known since early 80s, when the spectrum feature, ankle., was discov- 
ered in the Haverah Park data [23] at £^ ~ 1 x 10^^ eV. This feature was 
interpreted as transition from galactic to extragalactic cosmic rays (in con- 
trast to the calculations above where the ankle naturally appears as a part 
of the dip). The hypothesis of the transition at the ankle can be described 
phenomenologically as follows: At energy below 1 x 10^^ eV the cosmic ray 
flux is galactic and above - extragalactic. The galactic spectrum can be taken 
basically as power-law oc E~^^'^^ , but agreement with observations needs steep- 
ening at E > Eg, described by some steepening parameter. In effect one can 
use the parametrization: 

7gai(£;) = isTgai^"^^^' [l - a exp(6 log E/W'^ eV)l . (3) 
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The extragalactic generation spectrum is assumed to be power-law with index 
7g. Together with two constants of the normahzation for both fluxes, one has 
as minimum six free parameters to fit the observed spectrum. We found the 
best fit shown in Fig. 6. It is characterized by — 9.1 for 19 energy bins and 6 
free parameters, i.e. by ^ jd.o.j. — 0.7 for d.o.f.=13. The value ^ jd.o.j. < 1 
for the Poisson statistics signals for the large number of the free parameters 
and for very good formal fit to the experimental data. The problem of this 
ad hoc model is whether there is a physical model for propagation of galactic 
cosmic rays, which results in spectrum given by Eq.(3). It is just assumed 
that the spectrum at < 1 x 10^^ eV is the same as observed, while the 
dip model predicts this spectrum in excellent agreement with observations. 
An intermediate case is given in [24] where the dip is mostly described by 
extragalactic protons interacting with CMB with a small correction given by 
galactic cosmic rays only at the low-energy, E Ki 1 x 10^^ eV, part of the dip 
(see Fig. 13 in [24]). 

How does extragalactic magnetic field affect the discussed spectra features? 
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Fig. 6. The ad hoc model for explanation of the dip. The distorted power-law 
spectrum (3) is shown by curve "galactic" . The extragalactic spectrum with 7g = 2 

is normalized by the data in the interval (4 — 6) x lO"*^^ cV. The parameters of the 
galactic spectrum are found to provide the best fit to the data by the total spectrum 
(galactic -|- extragalactic), shown by the thick curve. It gives 7gai = 3.28. 



The infiuence of magnetic field on spectrum depends on the separation of the 
sources d. There is a statement which has a status of the theorem [25] : 

For uniform distribution of sources with separation much less than charac- 
teristic lengths of propagation, such as energy attenuation length /att and the 
diffusion length /dis, the diffuse spectrum of UHECR has an universal (stan- 
dard) form independent of mode of propagation. 
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For the realistic intergalactic magnetic fields the spectrum is universal in en- 
ergy interval 1 x 10^^—8 x 10^^ eV [18]. Note, however, that generation spectrum 
is defined in [25] as one outside the source. In this work we implicitly assume 
that the sources are transparent for UHE protons and thus the generation 
spectrum is the same as the acceleration spectrum. 

The most probable astrophysical sources of UHECR are AGN. They can ac- 
celerate particles to E 

max 1 X 10^^ eV and provide the needed emissivity 
of UHECR Co^Sx 10^^ erg/Mpc^r- The correlation of UHE particles with 
directions to special type of AGN, Bl Lacs, is found in analysis of work [26]. 
AGN as UHECR sources in case of quasi (rectilinear) propagation of protons 
explain most naturally the small-scale anisotropy [27]. 

The UHECR from AGN have a problem with superGZK particles with ener- 
gies ii^ > 1 X 10^° eV: (i) another component is needed for explanation of the 
AGASA excess, and (ii) no sources are observed in AGASA and other arrays 
in direction of superGZK particles. These problems probably imply the new 
physics, such as UHECR from superheavy dark matter, new signal carrier, 
like e.g. new light stable hadron, strongly interacting neutrino, and Lorentz 
invariance violation. For the last case it is interesting to note that if Lorentz 
invariance is weakly broken for e"'"e~ production, but strongly for pion pro- 
duction, then the modification factor is given by the curve r]ee in Fig. 3. This 
prediction agrees well with the Akeno- AGASA spectrum. 
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